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AMENDMENTS TO THE SPECIFICATION 

Page 1, line 18, to Page 2, line 7, replace with the following: 

-Many difficulties are encountered in such a magnet design problem, and many solutions 
have been proposed. The various proposed solutions exhibit one or more disadvantages. For 
example, the problem of volume minimization of electromagnetic coil systems has been considered 
previously by Kitamura et al M. Kitamura, S. Kakukawa. K. Mori, and T. Tominaka. "An 
optimal design technique for coil configurations in iron-shielded MRI magnets." IEEE Tran. Magn., 
vol. 30, no. 4, pp. 2352-2355, 1994, ("Kitamura"), and by Xu m et al, H. Xu, S. M. Conollv, G. C. 
Scott, and A. Macovski, "Homogeneous magnet design using linear programming," IEEE Trans. 
Masn.. vol. 36, no. .2, pp. 476-483, 2000 ("Xu") . Th e numb e rs in squar e brack e ts ar e r e ferenc e s to 
pubHsh e d pap e rs that ar e id e ntifi e d in App e ndix B ann e x e d h e r e with. The contents of e ach of thos e 
pubhcations ar e incorporat e d h e r e in by r e f e r e nc e . For simplicity, when discussion of a particular 
reference is needed or desirable, it may be identified by the name of the lead author or as, for 
e xampl e , Ref 1, which will b e understood to mean the Kitamura paper, or R e f 2, which will b e 
und e r s tood to m e an th e Xu pap e r, and so on . The reader is urged if needed to review R e f 1 and 
R e f 2 these references for a more complete discussion of the various known methods and their pros 
and cons and for a more complete understanding of the methods described herein. While the 
methods disclosed by R e fs. 1 Kitamura and 2 2Cu have benefits, they also have requirements that are 
undesirable. For example, the method of Rei^ Kitamura assumes unidirectional currents in the 
coils, while the method of R e f 2 Xu requires that the length-to-width ratio of the coils be specified. 
Therefore, neither approach determines, in general, the minimum volume solution, as is defined in 
this invention. It is also noted that a variety of procedures for optimizing electromagnets have been 
described that utilize criteria other than coil volume minimization; for a review of these the reader is 
referred to R e f 2 Xu . — 

Page 3, lines 11-21, replace with the following: 

" As with those of R e fs. 1 Kitamura and 2 2Cu, the method described here employs linear 
programming (LP) computation using a computer. However, a straightforward application of linear 
programming may not, by itself, be a practical approach for accurately calculating minimum volume 
coil arrangements. The reason for this is that the memory requirements and computational time for 
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a linear programming computation increase rapidly as one increases the number of numerical grid 
points used to represent the feasible volume. See the discussion in R e f. 2 Xu . To circumvent this 
problem, I first use linear programming with a sparse grid to obtain a coarse grained solution. This 
solution is then refined by numerically solving with a computer a set of nonlinear equations. The 
solving of these equations on a dense grid is feasible, because the memory requirements and 
computational time increase only linearly with the number of grid points. ~ 
Page 5, lines 3-10, replace with the following: 

~ I will first define in Section II a statement of the problem in appropriate mathematical 
terms; I will next provide in Section III the mathematical background for one skilled in this art to 
understand my new method; I will then describe in Section IV my method for minimizing coil 
volume including various optional steps for certain purposes; I will then describe in Section V my 
method for minimizing power; next; I will describe in Sections VI and VII some examples and 
results in comparison with other proposed solutions and some remarks; and finish in Section VIII 
with some summarizing comments. App e ndic e s A, B, and C follow Following Secfion VIII are 
Section IX, labeled Mathematical Proofs, Section X. labeled Tables L IL and IIL and finally the 
claims and Abstract. ~ 

Page 6, line 12, replace with the following: 
-- The magnetic field generated by the magnet is given by D. B. Montgomery, Solenoid Masnet 
Desisn, New York: Wiley, 1969 ("Montgomery") - 

Page 8, last two lines, replace with the following: 

~ The method for determining the current density J is founded on three principal 
mathematical theorems. Here I state these theorems, leaving their proofs for App e ndix A Section IX 
below . -- 

Page 9, lines 12-14, replace with the following: 

-- 1 note that this formulation has some similarities to, but is not identical to, ones presented 
by Kitamura a t al. [1] and by Xu e t al. [2] . For instance, J is assumed to be nonnegative in R e f 1 
Kitamura , while the condition of Eq. (12b) is not utiKzed is R e f. 2 Xu . — 
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Page 10, lines 7-15, replace with the following: 

-- This theorem shows how the minimization of V with equality constraints can be reduced 
to solving a system of nonHnear equations for the Aj. These equations are obtained (as is discussed 
below) by combining Eqs. (13) and (14). A somewhat related result has been previously derived by 
me for permanent magnet structures {€] , J. H. Jensen, "Optimization method for permanent-magnet 
structures," IEEE Trans. Magn., vol 35, no. 6, pp, 4465-4472, 1999 ("Jensen") . However, the 
equations and assumptions of that earlier development for a permanent magnet structure are not 
applicable to the electrical coil arrangements which can be designed using the methods and concepts 
of the present invention which are not disclosed in R e f [6] Jensen. A useful corollary of the present 
invention is that a solution with Xo = 0 must have Jo equal to the smallest current density magnitude 
that is consistent with the conditions of Eqs. (8) and (13). — 

Page 11, line 9, to Page 12, line 5, replace with the following: 

~ With the introduction of the numerical grid, the form of the problem given by Theorem 1 
becomes a standard linear programming problem as described in R e fs. 4 and 5 G. Hammerlin and 
K.-H. Hoffman, Numerical Mathematics. New York: Springer- Verlag, 1991 ("Hammerlin"):and W. 
H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannerv , Numerical Recipes in C: The Art 
of Scientific Computing. New York: Cambridge University Press, 1992 ("Press") . The functions J 
and U may be regarded as N-dimensional vectors and the optimization problem has 2N degrees of 
freedom. The total number of inequality constraints, corresponding to Eqs. (5), (6), and (12), is 4N 

+ 2Mh + 4Mf. 

A variety of algorithms have been proposed for solving linear programming problems, with 
the most popular being the simplex method {4^ Hammerlin, Press . All of these are usable in my 
invention, but I prefer to use a double precision version of a simplex method subroutine given in 
Ref 5 Press and have found it to be effective when N is a few hundred or less. For larger values of 
N, three problems arise. First, the method requires sufficient memory for a matrix with 
approximately 8N^ elements or about 64N^ bytes. On many computers, this limits N to no more 
than a few thousand. Second, if £ is small, as is the case when the homogeneous region is required 
to have a high degree of uniformity, then numerical rounding errors become important. In practice, 
I have found that this often restricts N to be less than one thousand. Finally, the average 
computational time scales roughly as N"^, which again can make large values of N impractical, 
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To be sure, these difficulties are to some extent algorithm dependent. In principle, both the 
memory and rounding problems could be largely circumvented, but at the price of a longer 
computational time. However, no linear programming algorithm has been proposed with a 
computational time that, on average, grows substantially more slowly with N than the simplex 
method {4} Hammerlin . Thus, a straightforward application of linear programming to the 
determination of the ideal current density is likely to be limited to sparse grids.— 

Page 14, lines 3-5, replace with the following: 

~ There are a number of known standard numerical methods for solving equations of the 
form of Eq. (25) (See R e f [5] Press ). All of these can be used in accordance with my invention. 
However, for this particular problem, I have found an iterative approach, similar to one described in 
R e f [6] Jensen , to be especially effective. In order to employ this method, Eq. (25) must be recast 
into a different form. Let me define an (M + 1) x (M + 1) matrix- 
Page 16, last line to Page 17, Hne 2, replace with the following: 

~ This follows fi-om inspection of Eq. (14). The current density is then found by repeating 
Step 3, as described above. (A more detailed discussion of the imposition of auxiliary conditions is 
given in R e f 6 Jensen .) ~ 

Page 18, line 13, to Page 19, Hne 6, replace with the following: 

- In order to illustrate the optimization method, I apply it to two different examples. First I 
compute the minimum volume coil arrangement for a magnet that is comparable to a design 
proposed in the seminal work Garr e tt [7] M. W. Garrett, "Thick cvlindrical coil systems for strong 
magnetic fields with field or gradient homogeneities of the 6th to 20th order," J. Appl. Phvs., vol. 
38, no. 6, pp. 2563-2586, 1967 ("Garrett") . In this paper, Garrett gives numerous coil arrangements 
that generate uniform fields. However, these arrangements are not optimized in any systematic 
manner. I choose one simple example and demonstrate how my method can lead to an arrangement 
with the same external geometry, current density magnitude, field strength and homogeneous 
volume, but with a significantly smaller coil volume. I also consider varying the current density 
magnitude and show how to minimize the power consumption for resistive coils. As a second 
example, I minimize the coil volume for a magnet similar to conventional actively shielded 
superconducting magnets used for whole-body MRI. 
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The arrangements considered are all symmetrical with respect to the z = 0 plane. As a 
consequence, the expansion coefficients, aj and bj, vanish automatically for even j. The constraints 
corresponding to these even order coefficients may then be excluded from the nonlinear system, and 
all the Xj with even j, except Xq must be zero. While this symmetry simplifies the calculation of the 
ideal current density, it should be emphasized that it is not a requirement for the application of my 
method. -- 

Page 19, line 8-16, replace with the following: 

- Figure 3a shows a three coil arrangement 12, 14, 16, taken from Fig. 3 and Table XVI of 
the work of Garrett f?}. The total length of the magnet is 2.8839 lo, the inner radius is 0.75 lo, and 
the outer radius is 1 .25 lo , where lo is an arbitrary length scale. The total volume of the coils is V = 

5.0813 lo^ The magnitude of the current density is Jo = 4.9086 Bo/(^ol<)), where Bo is the desired 
field strength, and the current flows in the same direction in all of the coils. The homogeneous 
volume with a 20 ppm peak-to-peak field variation is approximately an oblate spheroid with major 

and minor semiaxes of 0.42 lo and 0.36 lo. The coil arrangement is designed to cancel the 

harmonics of Eq. (16) for 2 < j < 8. -- 

Page 20, lines 5-12, replace with the following: 

" To obtain a design with rectangular coils, I follow step 4 of the optimization procedure and 
partition the domain D as indicated in Fig. 4a by the 13 rectangles shown. This leads to the coil 
arrangement shown in Fig. 4b 20-40. In order to obtain a high degree of accuracy, I used N = 70, 

000. The volume of the coils is V = 4.62 to^ which is only slightly larger than that of the ideal 

solution. The design of Fig. 4b consists of 13 coils. Each coil may be specified by giving the 

coordinates of the point lying nearest to the origin, {r{^, Zin), and the coordinates of the point lying 

farthest from the origin, (rout> Zout)- These values are listed in Table I of App e ndix C below for the 

coils that lie in the positive z half-plane. - 

Page 20, last line, to Page 21, line 6, replace with the following; 

- I now consider optimizing a magnet 110 with an inner radius of 50 cm, an outer radius of 
75 cm, and a total length of 125 cm as shown by the rectangle of Fig. 7b. I set Bo = 1.0 T and require 



{W:\05986\OOOH648000\003781 50.DOC liiiiililllilllliiJIililil } 



6 



Application No.: 09/732,322 Docket No.: 05986/OOOH648-USO 

the homogeneous volume to have a diameter of 45 cm, with a peak-to-peak field variation of less 
than 20 ppm. The magnitude of the fringe field 1 12 is restricted to be less than 0.5 mT (5 G) for 
distances greater than 400 cm from the center of magnet. These specifications are typical for a 
compact superconducting magnet designed for medical MRI, fS} S. Crazier and D. M. Doddrell 
"Compact MRI magnet design by stochastic optimization/' J. Magn. Reson., vol. 127. pp. 233-237, 
1997 r"Crazier"V 

Page 29, top line, replace with the following: 

- APPENDIX A Section IX-Mathematical Proofs - 

Page 30, 6^^ line from bottom, to Page 31, line 2, replace with the following: 

-- The function on the right side of Eq. (AlO) is analytic in the variable z. A standard theorem 

of analytic functions then indicates that if Eq. (AlO) holds for a dense set of z values, then it must 

hold for all z values {9] .R. V. Churchhill J. W. Brown, and R. F. Verhev, Complex Variables and 

Applications. New York: McGraw-Hill 1974. pp. 286-287 ("Churchill") . However, with the aid of 

Eqs. (3), (4), and (A2), one can show that the right side of Eq. (AlO) vanishes as z ^ oo. I then 

conclude that Eq. (AlO), and hence Eq. (A9), cannot be satisfied on a dense set of points, and 

therefore the volume associated with Db must be zero. From this it follows that = ^ (except 

possibly for a set points of zero measure) and that the condition of Eq. (12b) may be replaced with 

the stricter condition of Eq. (8). — 

Page 35, Delete the entire page. 

Page 36, top line, replace with the following: 

- APPENDIj ^ SECTION X 
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